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A review of supersymmetric dark matter in minimal supergravity unification 
' ■ ' I with R-parity invariance and with radiative breaking of the electro-weak sym- 

' metry is given. The analysis shows the lightest neutralino is the LSP over most 

of the parameter space of the supergravity model. The event rates in neutralino- 
lO ' nucleus scattering in dark matter detectors are also discussed. It is found that 

I the event rates are sensititive to the constraint from the 6 — > S7 experiment.lt 

(~| ' is also found that the event rates are sensitive to the constraints of relic density 

and in our anaysis we have used the accurate method for the computation of 
the neutralino relic density. Finally,the effect of the new results on quark polar- 
QJ \ izabilities, from the data of the Spin Muon Collaboration, on event rates is also 

• discussed. The analysis shows that the event rates for the Ge detectors and for 

other detectors which use heavy targets are only negligibly affected. 



■ 1. Introduction 



Considerable evidence for the presence of dark matter in the universe exists: in 
our galaxy, in other galaxies and in galactic clusters. The rotation curves of lumi- 
nous matter in spiral galaxies, point to massive amounts of non-luminous matter in 
the halo of galaxies and provide perhaps the strongest evidence for the existence of 
dark matter^. There are many possible candidates for dark matter both in particle 
physics and in astronomy. Thus in astronomy possible candidates for dark matter 
are Jupiters, brown dwarfs, neutron stars, black holes etc, while in particle physics 
one has the possibility of axions, neutrinos, sneutrinos, neutralinos etc. An important 
constraint arises from the fact that not all of the dark matter in the Universe can 
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be baryonic in nature. First, the baryonic dark matter is constrained severely from 
analysis of nucleosynthesis which show Qb < 0.1. Second, the recent results of the 
MACHO Collaboration^, and from EROS^ indicate that at best only 20 - 30% of 
the halo of galaxies is composed of MACHO'S (Massive Compact Halo Objects), and 
thus the remainder must be non-baryonic dark matter (NBDM). The non-baryonic 
dark matter could be either hot (HDM) or cold (CDM). The HDM could be one of 
the neutrino species (most likely possibility is the tau neutrino Ut), while the CDM 
could be either an axion, or a SUSY particle'' (a sncutrino i> or a neutralino x)- 
the following we shall pursue the SUSY possibility and assume that CDM is either 
a u or a X- Actually it turns out that in supergravity unification^, with radiative 
breaking of the electro-weak symmetry, the lightest neutralino turns out to be the 
LSP over most of the parameter space of the model^. Thus the model actually pre- 
dicts the neutralino to be the CDM^. Further, we assume a mix of cold and hot 
dark matter in the ratio of 2:1 as indicated by the COBE data. The quantity that 
appears in theoretical analyses is Qh"^, where h is the Hubble parameter in units of 
lOOKm/sMpc. Currently the experimental uncertainty in h is given by^ 0.82 ± 0.17 
(Freedman et. al) ; 0.53 ± 0.05 (Sandage et. al). For the purpose of our analyses 
here we shall assume an h in the range OA < h < 0.8 consistent with the above data. 
Then assuming that Q,b — 0.1, and ^cdm '■ ^HDM—'^'i one finds that 

0.1 < flh^ < 0.4 (1) 

As mentioned above there are two neutral states, the lightest neutralino x 
the sneutrino u, which are possible candidates for CDM in supersymmetric models, 
for example, the MSSM. However, in the MSSM there are many arbitrary parameters 
and the model does not predict what the LSP is. The situation is radically different 
in supergravity unified models^. Here the parameter space is five dimensional and 
reduces to a four dimensional space after fixing the Z-mass using radiative breaking of 
the electro-weak symmetry. The parameter space of the model is then characterized 
by mQ,mi/2, AQ,tan(3 where mo is the universal scalar mass, mi/2 is the universal 
gaugino mass, Aq is the universal trilinear coupling at the GUT Scale and tan/3 = 
V2/V1 where V2 gives mass to the top quark and Vi gives mass to the bottom quark. 
We limit the parameter space by the fine tuning criterion mo, mgiuino < ITeV where 
fngiuino — {otz/ aG)mi/2 and also limit tan(3 (since tan(3 is the ratio of two VeVs, a 
large tan(5 also implies a finetuning) so that tanP < 20. In this domain, we find that 
the lightest neutralino is also the LSP over most of the allowed region. 

2. Neutralino Relic Density 

The neutralino relic density at current temperatures is given by^ 

VLho^ = 4.75 X 10-^°{TzjT^f{T^/2.75fNl/\GeV~yj{xf))g/cm^ (2) 



Here is the current photon temperature,nir is the effective degrees of freedom 
computed at the freeze-out temperature, {T^/T^Y is a reheating factor, and J{xf) is 
given by 

J{xf) — J < av > dx (3) 

< av > is the thermal average of the neutrahno annihilation cross section and x = 
kT/rUz^. Now J receives contributions from Z-exchange, Higgs exchange and from 
the s-fermion exchange in the t-channel. For the s-fermion exchange the conventional 
approximation^ of expanding < av >= a + hu'^ in the integrand in eq. (3) is valid 
and we use this approximation. However, for the Z and Higgs pole terms such an 
expansion is a poor approximation in the region below the poles^"^^. Thus for Jmggs 
and Jz-poie we use a rigorous thermal averaging over the poles. For example, consider 
the annihilation via the hghtest Higgs pole. Here 

(^v = "^"'^^^ ^i— — — (4) 

= {ml - 4m^2/m^2 (5) 
-fh = 'mhTh/mz2 (6) 

where irih is the Higgs mass, and Th is the Higgs width. 

Computationally eq. (3) implies a double integration over a pole. Since the pole 
is associated with a very small width numerical integrations are tricky as a sharp 
resonance can be easily missed. A more reliable procedure is to reduce eq. (3) to a 
single integral so that-'^'^"^-'^ 

^H....(^/) = ^^|^[/i.+ |/2.] (7) 
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A similar analysis can be carried out for Jz- What one finds then is that eq. (7) and 
the similar expression for Jz give a smooth result on integration over the poles. 

3. Analytic Analysis of the Neutralino Composition 

The lightest neutralino Zi is a linear combination of the four neutral states 
{W,B,Hi,H2) so that 



Zi = niW^ + n2B + mH^ + 714^2° 



(10) 



where = 1—4) are the co-efficient of the components of the Zx eigenvector and are 
discussed below. The co-efficients play an important role both in the relic density 
analyses as well as in the analyses of event rates in ncutralino-nucleus scattering. It 
is useful then to gain an analytic understanding of the parametric dependence of ni 
on the basic paramters of the model. The rzj are determined by the neutralino mass 
matrix 
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where TTij = {ai/a^jnig, a — Mzcos9wcos(3, b — —Mzcos9wsin(3, c — —Mzsin9wcos(5 
and d = Mzsin6wsin(5. Here ifi2{nii) are the SU(2)(U(1)) gaugino masses deter- 
mined by the relation fhi = {ai/ aG)f^i/2- The parameter /i is determined by fixing 
the Z-mass using radiative breaking of the electro-weak symmetry. It is found that 
over most of the parameter space n is determined to be large i.e. |/i^/M|| >> 1. In 
this domain one can carry out a perturbative expansion of Ui in Mz/i^- One finds 
that to second order in [Mz/ jj) one has-*^^ 
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The expansion of eqs. (12-15) is found to be accurate to (3 — 5)% over a signifi- 
cant region of the parameter space. From the above one can easily see that n2 > 
771,773,774, \n-i\ > \n4\ whcrc in getting this result we have used the radiative electro- 
weak symmetry breaking relation tanf3 > 1. We note that eq. (13) implies that the 
neutralino is mostly a Bino in the scaling limit. However, a note of caution is needed 
in that one should not take the \Mz/n\ oo hmit. This limit is dangerous since the 
gaugino- higgsino interference terms which are proportional to (77i, 772). (773, 714) vanish 
in this limit. In realistic analyses, as in the computation of the coherent part of 
neutralino-nucleus scattering, such terms make significant contributions and cannot 
be set to zero. 



4. Neutralino Detection via Neutralino - Nucleus Scattering 



Various possibilities for the detection of neutralino dark matter have been dis- 



cussed in the literature. For example, annihilation of neutralinos in the galactic halos 
can produce an observable signal, i.e., % + x — A + X, where A can be an energetic 
gamma ray, positron or an anti-proton'^. However, the backgrounds in these processes 
are rather significant so this process does not seem very encouraging for the detection 
of the neutralino. A second possibility is that the halo neutralinos are captured in 
the center of earth and sun, annihilate and produce upward moving neutrinos (and 
muons) in detectors on earth^^'^^'^^. The background in these processes are signifi- 
cantly reduced due to the angular windows around earth and sun. Current estimates, 
however, show that one needs around 10^ — lO'^m^ neutrino telescopes to see any 
significant effect. The telescopes currently being planned aim to approach O(lO^m^) 
area. So once again this mode of detection also does not appear very optimistic for 
the neutralino dark matter. The most optimistic mode for detection of neutralinos 
appears to be scattering of neutralinos off nuclei in cryogenic detectors^'^'^^'^^'^^'^^'^^. 
Several detectors using this mode are in various stages of development. We shall focus 
on this mode of detection for the rest of the talk. More recently there has also been 
a discussion of detection of neutralinos via atomic excitations^^, but at the moment 
this possibility requires more investigation. The prime detector in neutralino-nucleus 
scattering is the quark and the effective interaction that governs the neutralino-quark 
scattering consists of a spin-dependent(incoherent) part and a spin-independent (co- 
herent) part, and is given by^ 



Ceff = Z I'y ^-^5 ZiqY{AqPL + BqPji)q + ZiZirriqqCqq (16) 

Here Aq, Bq are the spin-dependent amplitudes which arise from the s-channel Z- 
exchange and the t-channel squark exchange, and Cq is the spin-independent am- 
plitude arising from the s-channel Higgs exchange and t-channel squark- exchange. 
Realistically, of course the quarks are bound inside nuclei so a reasonable amount of 
nuclear physics enters in the anlaysis. The total event rate is given by^ 



R — [Rcoh + Ry 
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(18) 



(19) 



In the above p^i is the local density of dark matter v^.^ is the relative velocity, 

— * 

is the nucleus mass, J is the nucleus spin and A is defined by A < A^|J|A^ >=< 

— * 

A^|ES'j|A^ > where S-i is the nucleon spin. From eqs. (18) and (19) one finds that 
for large Mjv one has Rcoh ^ 0{Mn), Ri^c ~ 0{X^J{J + 1)/Mn). Thus in principle 
one has two qualitatively different types of detectors, i.e., those with large Mn and 



those with large values of J( J+ 1). In practice there is seldom a case where Rcoh is 
negligible and realistic analyses even for light target material (e.g.,^He,CaF2) require 
inclusion of both i?j„c and Rcoh- For heavy targets (e.g., Ge,NaI,Pb) Rinc is typically 
small, i.e. only a few percent of the total R. 

We note that there are uncertainties in the evaluation of both R^nc and Rcoh- 
Uncertainties in R^nc arise due to experimental uncertainties in the determination 
of Aq on which Rinc depends, where Aq are the quark polarizabilities defined by 
< p(n) 1 57^755 |p(n) >= S'jf'^'^Aq, where 5*^*^"^ = (0,5^'^")) is the nuclcon spin. We 
shall discuss later the sensitivity of the results to the determination of Aq. There is 
also an uncertainty in the determination of Rcoh- This arises due to the uncertainty 
in the determination of s-quark matrix elements^^ < n\msSs\n >= fsM^ that enter in 
the computation of C^. Currently fg has a significant uncertainty, about 50%, which 
can lead to an uncertainty of 0(30 — 50)% in the determination of Rcoh- 

Next we discuss briefly the relative contribution of the heavy neutral Higgs to the 
coherent part of the scattering. Naively one might expect the heavy Higgs contribu- 
tion to be neglible since it would be suppressed by the heavy Higgs {mass)^ while 
the light Higgs contribution is suppressed only by the light Higgs {massy. However, 
this assessment is correct for the up quark but not for the down as can be seen by 
the expression for Cg below: 



Here a is the Higgs mixing angle. Now the Higgs mixing angle a is typically small 
so from eq. (20) one finds that there is a suppression of the d-quark contribution in 
the light Higgs sector which often can be more than the {massy suppression in the 
heavy Higgs sector. Thus the heavy Higgs contribution cannot be neglected as it can 
make a substantial contribution to Cg. 

5. b ^ s^ Branching Ratio Constraint on Event Rates 

Recently the CLEO Collaboration obtained the first experimental determination 
of the photonic penguin process b — > 57. For the inclusive b ^ sj decay they find 
the result^^ 




(20) 



where 



Fh = {ni — n2tan9w) {sinans + cosan^) 
Fh = {ni — n2tan9w){—cosans + sinan^) 



(21) 
(22) 



BR{b 57) = (2.32 ± 0.51 ± 0.29 ± 0.32) x 10 
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(23) 



Now in the SM 6 — > 57 decay receives contributions from the W-exchange. For the 
super symmetric case there are additional contributions arising from the exchange 
of ,W , Z ,and g. There are many uncertainties in the theoretical evaluation of 
b 57. These uncertainties arise from uncertainties in the determination of ag, 
quark masses KM matrix elements, and uncertainties arising from the next to leading 
order (NLO) QCD corrections which can be as large as 0(30)%. There is a significant 
debate in the literature currently regarding what exactly the BR{b — > 57) value is in 
the SM^^. An accurate answer to this question can only result after the next-to-lcading 
order(NLO) QCD corrections have been computed reliably. Similar uncertainties are 
present in the evaluation of 6 ^ 57 in SUSY theory. The discussion of 6 ^ 57 
constraint on neutralino relic density ^^'^^ can be facilitated by use of the parameter 
rsusY^^ which we define by the ratio tsusy — BR{b — > s^)susy / BR{b — > S'j)sm ■ 
Many of the uncertainties discussed above cancel out in the ratio tsusy- However, 
we must keep in mind that the NLO corrections would in general be different for 
the SUSY case than for the SM case. In this analysis, however, we limit ourselves 
to the leading order evaluation. Analogous to tsusy we can also define^^ T^xpt = 
BR{b — > s'y)expt/BR{b — > s'y)sM ■ The CLEO Collaboration use an SM value of 
(2.75 ± 0.8) X 10-^. Using this value and the result of eq. (23) one finds that Vexp lies 
in the range rexp = 0.46 — 2.2. This range turns out to be a rather strong constraint 
on SUSY theory if we assume that tsusy — fexp^^- This is so because in the SUSY 
case one normally gets a much larger range for tsusy , i-e., a range of ~ (0, 10). 

6. Analysis and Results 

We first discuss the event rates without inclusion of the b ^ s'j constraint. Results 
on maximum and minimum of event rates are exhibited in Fig. 1 as a function of 
T^giuino for the target materials: Pb, Ge and CaF2. We note that the maximum curves 
of the event rates all exhibit a dip when rugiuino is in the mass range ~ 300 — 450GeV 
independent of the target material. This dip arises^^ due to the Z-pole and the Higgs 
pole effects and the relic density constraint. Effectively, rapid annihilation near the 
Z and Higgs poles leads to values of Vth"^ which fall below the CHDM limit and have 
to be eliminated. The eliminated part of the parameter space contains the light 
SUSY spectrum and the large event rates. Thus one sees a dip in the event rates in 
the vicinity of m^^ ~ Mz/'2, and ~ irih/'^- The total event rate R is generally 
dominated by Rcoh for all except the lightest target materials. Now Rcoh depends 
on the gaugino-higgsino interference term which is proportional to (ni, 712) x (^3, 714). 
Using the behavior of rij with large ji and the dependence of Rcoh on the gaugino- 
higgsino interference, we can understand the behavior of R for large irigiuino- Typically 
/i is an increasing function of rrigiuino, and the n,(i=l,3,4) are decreasing functions as 
H increases. Thus Rcoh falls monotonically^^ with nigiuino beyond the dip as can be 
seen in Figs. 1. (A similar analysis holds as a function of mo). Now it is easily seen 



Fig. 1. Plot of maximum and minimum event rates as functions of nigiuino for CaF2 (dash-dot), 
Gc(dash) and Pb(solid) when tanfS < 20 and other parameters span the parameter space.. Fig. la 
is for /X < and Fig. lb for /z > 0. 



Fig. 2. Same as Fig 1 except that the b —^ sj constraint of rsuSY < 1-33 is imposed. 

that for nigiuino > 650GeV for < and nigiuino > TQQGeV for /x > the event rates 
even for heavy targets (e.g. Pb) fall below the level of 0.01 event/kg. d. This is the 
level of sensitivity that detectors in the next 5-10 yrs hope to achieve. 

Next we discuss the effect of the 6 — > 57 constraint on dark matter^^'^^'-^^. The 
quantitivc effects of 6 — > 57 depend on the value of tsusy- Generally, the b ^ 
constraints is more severe for fi > than for fi < 0. Similarly, one finds a significant 
effect on the event rates for fi > 0, while the effect for < is relatively smaller. 
Results for the case tsusy < 1-33 is shown in Fig. 2 for the cases > (Fig. 2a) 
and jJL < 0(Fig. 2b). As indicated above the allowed region of the parameter space 
shrinks significantly for > (Fig. 2a) and the maximum allowed event rates also 
fall. The corresponding effect on < (Fig. 2b) is signficantly smaller. 

Finally we discuss the effect of the variations in quark polarizabilities Ag on 
the event rates. The part sensitive to Ag is Rinc- The previous determinations of 
Ag using the EMC data gave^^, (Am, Ad, As) = (0.77 ± 0.08, -0.49 ± 0.08, -0.15 ± 
0.08). Recently there has been a reanalysis of Aq using new data from SMC^^ which 
gives^^ {Au,Ad,As) = (0.83±0.03, -0.43±0.03, -0.10i0.03). These determinations 
are consistent with each other within la, but the variations of Aq, specifically the 
variation of As, can generate signficant changes in Rinc- However, i?j„c is generally 



a small component of the total R^^'^^. We have analysed several target materials, 
Hc,CaF2,Ge,Pb etc., and find that except for the lightest target materials, i.e. He 
and CaF2, Rinc is only a few precent of the the total. Thus the event rates for targets 
such as Ge are not significantly affected by the new determination of Aq^^'^^. 

7. Concluding Remarks 

We have given here an analysis of neutralino dark matter within minimal super- 
gravity unification. Remarkably the model predicts that the lightest neutralino is also 
the LSP over most of the parameter space of the model, and thus the model gives 
a candidate for cold dark matter. We have analysed the relic density of neutralinos, 
using the accurate method which integrates over the Z and Higgs poles in thermal 
averaging of av. These effects are found to be significant for values of gluino masses 
< 400^6^. We have analysed the event rates in neutralino-nucleus scattering and 
find that there is a significant region of the parameter space where event rates > 0.01 
event/kg.d are predicted for targets such as Ge and Pb. This region of the parameter 
space will be accessible to current and future technologies over the next 5-10 years. 
However, more sensitive detectors, two to three orders of magnitude more sensitive, are 
needed to probe most of the parameter space of the minimal supergravity model. 
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